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We study the evolution of a photonic state in a coupled pair of anharmonic cavities and compare 
it with the corresponding system of two coupled harmonic cavities. The photons in the anharmonic 
cavities interact with a two-level atom and are described within the Jaynes-Cummings model. Start- 
ing from an eigenstate of the Jaynes-Cummings model with A'^ photons, the evolution of this state in 
the presence of photon tunneling between the cavities is studied. We evaluate the spectral density, 

- - - the dynamics of the return as well as the transition amplitudes and the probability for the dynamical 

Cn ■ creation of a NOON state. 
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Ph, I. INTRODUCTION 

< 

Q>^ ' The dynamics of isolated many-body quantum systems has been a subject of intense research in atomic 

physics during recent years, in experiment [IJ as well as in theory 0. This interest has two major aspects. 
One is to prepare the system in a well-defined initial state and, secondly, to study its evolution for a 
period of time due to quantum tunneling and particle-particle interaction. The preparation of the initial 
state as the groundstate of a certain Hamiltonian Hq and the evolution with exp(— iift) for a different 
2 ' Hamiltonian H involves a sudden change Hq -^ H, which is usually called a quench. Such a quench can 

C^ , be realized in atomic systems by changing the potential wells in which the atoms are trapped [3|. For 

^ ■ instance, bosonic atoms are prepared in a Fock state, where a definite number of atoms are localized in 

O^l deep optical potential wells. Then the potential barrier between a pair of neighboring wells is suddenly 

reduced such that the atoms can tunnel between these wells [J, |5| . From the theoretical point of view the 
statistical properties of this problem have been studied intensively using the Hubbard model and related 
^ . models |6l4l3| . An essential element of the dynamical analysis is the evaluation of the spectral weights 

•/^ ' with respect to the initial state [ijj UM ' which links directly the quantum evolution with the spectral 

"^ , properties of the underlying Hamiltonian. 

Besides systems of ultracold atoms an alternative approach of controllable bosons is based on photons. 
Then the role of the potential wells is played by microwave or optical cavities. The experimental prepa- 
'nT ■ ration of Fock states in a optical cavity has been achieved recently [Ig, [iJ] ■ This is a crucial step towards 

^^ ' a systematic study of correlated many-body systems with photonic states. The interaction between the 

photons is indirectly mediated by atoms inside the optical cavities, which interact directly with the pho- 
tons |18l - l20l | . Once a Fock state with N photons has been prepared inside an optical cavity, we can couple 
the latter with another optical cavity by a waveguide or an optical fiber. Then the photons can tunnel 
between the two optical cavities, leading to a quantum evolution of the initial Fock state |-/V, 0) within 
the Hilbert space that is spanned by the eigenstates of the Hamiltonian of the new system. This new 
C^ ■ Hamiltonian can be approximated, for instance, by the Hubbard Hamiltonian, as suggested recently by 

several groups |2ll424l | . This type of system, including atomic degrees of freedom, was studied within a 
Hartree-Fock approximation [25|. 

The evolution from the initial Fock state can, in principle, lead to an entangled state, such as the 
NOON state {\N,0) + \0,N))/V2. The NOON state has attracted much attention because it can be 
used for highly accurate interferometry and other precision measurements |26l - l29j and for technological 
application such as optical lithograph y f30|. Various methods for the creation of photonic NOON states 
have been suggested in the literature [3l|, Is^ and indeed experimentally created up to A^ = 5 photons 
|29| . This clearly indicates that the creation of these entangled states is realistic. However, it still remains 
a problem to created NOON states for large N. The author discussed recently the dynamical creation 
of a NOON state from ultracold bosonic atoms in a double well [33] , which shows that a balanced effect 
of inter-well tunneling and intra-well interaction can produce such a state with moderate probability. 
Since photon-photon interaction can also be mediated in a cavity by coupling the photons to atoms, we 
will analyze in the following the dynamical creation of entan gled (NOON) states for a pair of anharmonic 
cavities, described by coupled Jaynes-Cummings models |34ll35|. 

The paper is organized as follows: In Sect. [H] we introduce the basic quantities for the dynamics of 
our system with two cavities. Then we discuss the return probability, the transition probability and their 
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relation with entangled states in Sect. IIIIl applying it to a single cavity with a two-level atom (Sect. 
IIII Al) . to a coupled pair of harmonic cavities (Sect. IIIIBp and eventually to a coupled pair of anharmonic 
cavities (Sect. IIIICp . The results of the recursive projection method of Sect. IIII CI are discussed and 
compared with the results of the coupled harmonic cavities in Sect. IIVI Finally, we summarize the work 
in Sect. E 

II. SPECTRAL DENSITY AND THE EVOLUTION OF ISOLATED SYSTEMS 

We consider a system which is isolated from the environment. In terms of photonic states this can be 
realized by an ideal optical cavity. With the initial state |^o) we obtain for the time evolution |^t) = 
e~'^^'^\'^o) or the evolution of the return probability K'l'ol^t)!^ with the return amplitude (^ol^t) = 
(^o|e~*^*|^o)- A Laplace transformation relates the return amplitude with the resolvent through the 
identity 

(*o|**> = J {^o\{z - H)-'\^o)e-'''dz , (1) 

where the contour F encloses all the eigenvalues Ej (j = 0, 1, ..., N) of H, assuming that the underlying 
Hilbert space is iV + 1 dimensional. With the corresponding eigenstates \Ej) the spectral representation 
of the resolvent is a rational function of z: 

{^o\iz-H)-^\^o)=f2^-^^^^J^ Q.^,{z) = l[iz-E,), (2) 

]-^ z-^j Qn+i{z) ^J-^ 

where Ppf{z), Qpf^i{z) are polynomials in z of order N, N + \, respectively. These polynomials can 
be evaluated by the recursive projection method (REM) [36|. The method is based on a systematic 
expansion of the resolvent (^oK-z ~ ^)~^|*o): starting from the initial state |^o)- It can be understood 
as a directed random walk in Hilbert space, where each subspace is only visited once. The latter is the 
main advantage of the REM that allows us to calculate efficiently the resolvent (\E'o|(-z — H)^^\'^q) on an 
N + 1-dimensional Hilbert space. 

The expression in Eq. ([2|) suggests the introduction of the photonic spectral density pofi{E) as the 
imaginary part of the resolvent: 

PoAE) = \im-Im{^o\iE~^e~H)-'\^o) - 1™ " E 2 | r F '^L, =T.\('^o\Ej)\H{E~ E,) , (3) 

j=Q "• 1' j=0 

where \^q) is a reference state. In other words, pq,o{E) is the diagonal element of the density matrix 
with respect to |^o)- Then the return amplitude can be written as the Fourier transform of the photonic 
spectral density 

N 

(*o|^t) = / po,o{E)e-''''dE = Y. \{^o\E,)\\-''^^' . (4) 

We can also evaluate other elements of the density matrix, such as the off-diagonal element 

Pi,o{E) ^ - lim Im{^>i\iE ^ le - H)-'\^o) =y^{'i'i\E,){E,\^o)S{E - E,) , (5) 

J 

whose Fourier transforms gives the transition amplitude between the states |^o) and |^i): 

N 

(*i|*t> = / Pi,oiE)e-''''dE = J2('^^\^o){Ej\^o)e-'''^' . (6) 



j=o 



To characterize the entangled state |^t) — cqI^o) +ci |^i) that may appear during the evolution, we need 
to evaluate the amplitudes (0]), ^ and count how often they realize certain values cq, ci simultaneously 
during a long period of time. After normalization, this defines the conditional probability P(co,ci) for 
having (^ol^t) = co and (^ij^I't) = ci at a given time t. 



III. RETURN PROBABILITY, TRANSITION PROBABILITY AND ENTANGLEMENT 

The central idea is to prepare an eigenstate of the cavity as initial state and then change the conditions 
of the system, either by adding a two-level atom to the cavity or by coupling another cavity through an 
optical fiber. As a result of this change, the system starts to evolve in Hilbert space to visit all possible 
eigenstates of the new Hamiltonian which have a non- vanishing overlap with the initial state. During its 
evolution the system may visit entangled states with certain probability. The dynamics and the entangled 
states will be calculated in the following subsections for three different cases. 

A. Single cavity with a tw^o-level atom 

An anharmonicity in a cavity can be created by adding an atom which interacts with the photons 
[ISl [l9l |37| . In the case of a single two-level atom we can describe the absorption and emission of photons 
by the atom approximately with the Jaynes-Cummings model [34, 35], whose Hamiltonian reads 

Hjc — woa^a + (wq + A)ctc - g{a)c + c^a) . (7) 

A is the detuning between the atomic excitation energy and the photon energy, c' (c) is the creation 
(annihilation) operator of the atomic excitation, and q is the coupling strength between the photons and 
the atom. The eigenvalues of this Hamiltonian are [3J, l35l | 



K.± = wo(^ + 1/2) ± v/A2 + 4g2(n + 1) (8) 

with eigenstates for A = 

|n,±) = ^=[±|n : 1) + |n + 1 : 0)], E^,± = wo(n + 1/2) ± 2.g\Am , 

where \n : j) is a Fock state with n photons and an atomic state with j — Q (atomic groundstate) and 
j = 1 (excited atom). Thus, the eigenstates are superpositions of two Fock states, one with n photons and 
the atom in the ground state and one with n — 1 photons and the atom in the excited state. This implies 
that the energy levels can be doubly degenerate with respect to the number of photons. A superposition 
of many eigenstates states for the initial state \^o) can lead to a more complex behavior, such as a 
collapse and revival dynamics [3a, [33] . 

The eigenstates of a harmonic cavity (without the atom) are Fock states \n) with n photons. We 
prepare a harmonic cavity in one of these eigenstates. This can be achieved in a real system, as recent 
experiments have demonstrated [la, [l3- Then we add a two- level atom in the ground state ]0), such 
that the initial Fock state of the combined system is a product state \n : 0) = 1«)10), whose evolution is 
described within the JC model as 

]*,) = e-'"-'^'\n : 0) = er'^'-^-^' ^\n - 1, +) + e-^^-^-*^]?! - 1, -) . 

Here we have assumed that the two-level system and the cavity are in resonance (i.e. A = 0) in order 
to have an optimal exchange between the atom and the photons. This implies for the return (transition) 
amplitudes 

(n : Ol*t) = e-''^o{n-i/2)t (.os{2gy/^t), (n - 1 : O]**) = _je-'"°("-i/2)« sin(2gV^i) 

and for the return (transition) probabilities 

\{n : Q\^t)? = cos2(2gV^i), \{n - I : l\^t)? - s\YV^{2g^t) . 

Thus, we see Rabi oscillations between the two Fock states \n : 0) and ]n — 1 : 1) with frequency 
ilfl — 2g^/n. This is an extreme case of Hilbert-space localization, where the system is constraint to 
a two-dimensional subspace. It is enforced by the fact that the eigenstates of the JC model are linear 
combinations of only two Fock states. The drawback of the extreme localization in Hilbert space is that 
we are not able to create dynamically entangled Fock states, except for the superposition of \n : 0) and 
]n — 1 : 1). In the general case the eigenstates of the Hamiltonian may be a superposition of many Fock 
states. Then the overlap of the eigenfunctions with the initial Fock state plays a crucial role. As a simple 
example we consider in the next section two harmonic cavities which are coupled by an optical fiber. 



B. Two coupled harmonic cavities 

The Hamiltonian of two uncoupled harmonic cavities is Hhc = <^o ^n=i 2 'A^j' where the index j = 1, 2 
refers to the two cavities, has product Fock states \N — k,k) = \N — k)\k) {k = 0, .., N) as eigenstates. 
Now we couple these cavities and obtain the Hamiltonian 

Hhc = -J{a\a2 + alai) + uJo{alai + 0302) 

with eigenstates |A^ — k;k). Here we assume that N is even and calculate the overlap of the eigenstates 
\N ~k;k) with the initial Fock state |iV, 0) [20l |: 

{N,0\N-k;k}=2-^/'i^^j , (9) 

which is non-zero for all eigenstates. The density-matrix elements with respect to |^o) — l^^O) and 
1*1) = \0,N) are 

^ /N\ ^ /N\ 

Po,o{E) = 2-^ ^ J<5(i? + J{2k - N), pME) = 2"^ ^ J {-1)'S{E + J{2k - N)) . (10) 

fc=0 ^ ^ k=0 ^ ^ 

Thus, there is a binomial distribution for the spectral weight |(iV, 0|A^ — k; fc)P, with a maximal overlap 
for an equally distributed number of photons. For large N the binomial distribution becomes a Gaussian 
distribution, where the width of the envelope is related to the energy level spacings A£' = 2J. The 
Gaussian result resembles the Central Limit Theorem for independent photons. Such a behavior was also 
found previously for freely expanding bosons from an initial Fock state (40[. A Fourier transformation 
reveals a periodic behavior of the return and transition amplitudes as 

(iV,0|e~'^*|iV,0) = cos^(Ji), (0,iV|e"*^*|iV,0) = (-z)^ sin^(Jt) . (11) 

Thus the evolution of the Fock state is periodic with period 2tt/ J but leads to a NOON state co|A^, 0) -|- 
cjv|0, A^) only with a probability that decays exponentially with N . For larger values of N the probability 
P{co,cisi) indicates an anti-correlation: P{cq,cn) vanishes as soon as both cq and cn become nonzero. 
Therefore, the overlap of \^t) with a NOON state is strongly suppressed. This is a consequence of the fact 
that for an increasing A'^ the particles disappear in the (A^ -I- l)-dimensional Hilbert space because there 
is no constraint due to interaction. 



C. Two coupled cavities with tw^o-level atoms 

Now we prepare an anharmonic cavity of Sect. IHI Al in the eigenstate of the JC model |A, +) and 
connect it with another anharmonic cavity which is in the state |0, -I-). After the connection the photons 
start to tunnel between the two cavities. This system is now described by the Hamiltonian 

H2JC = -J{a\a2 + 4ai) -t- ^ [a;oa]aj + woc]cj - g{a}jCj -I- c^a^)] , (12) 

J = l,2 

where the first term describes the tunneling of photons between the cavities with rate J and the second 
term represents the absorption and emission of photons by the two-level atom inside each cavity. For 
the initial state we prepare a product of JC eigenstates |A^, a;0,a-') = \N,a)\0,a'). The operators of the 
Hamiltonian (fT2|) act on the JC eigenstates separately. In particular, photon tunneling is controlled by 
the following matrix elements 

,, _, , ,, , VkTl + Vk ,, ^ , ,, ^ ^/k+l-^/k 
(fc- l,+|a|fc,+) = ^ , {k-l,+\a\k,~) = ^ , 

(fc-l,-|a|A:, +) = , (A;-l,-|a|fc,-) = 



(fc-l,+|at|fc,+) = 



Vk + 2 + Vk + l 


Vk + 2 


2 

-Vk + 1 



(fc-l,+|at|fc,-) = 



(fc-l,-|at|fc,+)= ^"^" ^"^\ (fc-l,-|at|fc,-)= v^ . -^ v^ , ^ ^ ^^3^ 



^/k + 2-^/k + l 


2 


Vfc + 2 + V^' + 1 



With the ratio 



Vk + 2-y/k + l 



y/k + 2 + y/k + l 



(14) 



the flipping of the atomic levels during the photon tunneling between the cavities is strongly suppressed. 
Moreover, we can use ^/k + 1 + y/k « 2^/k. With these approximations we decouple the ± states in the 
cavities to obtain a Hubbard-like model, where the n^ interaction is replaced by a ^/nj photon-photon 
interaction: 



H, 



eff 



- J(a{a2 + alai) + a;o(a{ai + 0302) + 2ag{ya[ai + y a\a2) ■ 



(15) 



Just like the Hubbard model, this Hamiltonian has a two-fold degeneracy for J = due to the equivalence 
of the two cavities. On the other hand, the interaction is weaker than the n^ interaction of the Hubbard 
model. This indicates that the properties of the coupled JC models may resemble the behavior of the 
Bose-Hubbard model in a double well [Sj], with less pronounced interaction features though. 

The appearance of i?e// brings us in the position to apply the RPM of Ref. [33] , only replacing the 
interaction term. Assuming that N is even, all projected spaces 7^2^ are two-dimensional and spanned by 
{|A^ — J, (t; j, cr), |j, cr; iV — j, cr)} (j = 0, ..., N/2). This leads to a recurrence relation in the base of the two 
JC states |A^, a\ 0, u), |0, a\ N, a) as initial states. The value of u = ±1 affects only the sign of coupling 
between cavity photons and the two-level system. Therefore, we ignore subsequently the a dependence 
in the matrix elements. If we define 



aN/2 ^ {N,0\{z - H)-^\N,0) = {0,N\{z - H)-^\0,N) 



and 



bN/2 = {0,N\iz-H)-'\N,0)^{N,0\iz~H)-'\0,N) , 
aN/2 and &7V/2 are obtained from the iteration of the recurrence relation (for details cf. |33|) 

1 fl 



_ , Ofc+i bk+i 
9k+i — \ 1 „ } , 9o — 7 \ n 



I] {k = 0,1,..., N/2 -I) 



with coefficients 
ak+i = 



z - fk+i - J^ak{N/2 + k + l)(N/2 - k) 



z - fk+i - J^akiN/2 + k + l){N/2 ^ k) - J%l{N/2 + k + l)2(iV/2 - k) 



(16) 
(17) 

(18) 
(19) 



bk+i 



and 



J^bk{N/2 + k + l){N/2-k) 



z - fk+i - .Pak{N/2 + k + l)(N/2 ^ k) - J^bl{N/2 + k + If {N/2 - k) 



fk+i = 2agy^N/2 + fc + 1 + 2agy^N/2 ~k~l . 
The recurrence relation terminates after N/2 steps with 

_ /ajv/2 6jv/2 
^""/'-[bj./^ a^i2 



(20) 



(21) 



(22) 
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Here it should be noticed that there exists an invariance of the recurrence relation under the following 
simultaneous sign changes in Eqs. P^ and (^01 

z — >■ — z, g — ^ — g, flj — )■ — Cj, hj — ^ — 6j . (23) 

This implies that a change from cr = + to cr = — in the initial JC states results in a mirror image with 
respect to energy of pQ^o{E,a) and pN,a{E,a): 

Po.aiE, a) = /9o,o(--E, -cr), PNfl{E, cr) = pNfl{-E, -a) . (24) 

Moreover, the density-matrix elements are invariant with respect to the harmonic frequency luq of the 
cavities, except for a global energy shift. This reflects an important universality of the density matrix 
that allows us to separate the harmonic from the anharmonic properties of the cavities. 

IV. RESULTS 

The properties of two coupled anharmonic cavities in Sect. IIII CI are characterized by two equivalent 
JC models and tunneling of photons between them. According to Eqs. p^ . pT)) . the iteration of Eqs. 
((T9| , (|20|) gives us the following four matrix elements of the resolvent 

{N,0\{z-H)-'\N,0), (0,iV|(z-i?)-i|0,iV), {0,N\{z - H)-'\N,0) = {N,0\iz - H)-'\0,N) . 

Moreover, according to Eq. ([2]) these matrix elements are rational functions of z. For N photons these 
are lengthy expressions with A^ + 1 poles. Therefore, it is convenient to present the results as plots with 
respect to the energy. 

Without inter-cavity tunneling the many-photon spectrum has a two-fold degeneracy due to the equiv- 
alence of the cavities. This degeneracy is lifted by the tunneling term, as one can see in the spectrum 
presented in Fig. [TJ Now we can compare this with the situation of two coupled harmonic cavities, as 
described in Sect. IIII Bl to evaluate the role of the photon-photon interaction. We start with the case 
of disconnected cavities (J = 0) and realize that the spectrum for N — k photons in one cavity and k 
photons in the other cavity is completely degenerate for harmonic cavities 

EN~kM = wo(A^ - fc) + ujQk = uj^N 

but only two-fold degenerate for anharmonic cavities 

EN-k,k = u}oN + ag{VN - k + \/k) . 

After connecting the cavities the degeneracy is completely lifted and an equidistant spectrum appears with 
level spacing Ai? — 2 J for the harmonic cavities in Eq. (jTOl) . J 7^ also lifts the two- fold degeneracies 
of the anharmonic cavities, as depicted in Fig. [T] However, the levels are more irregularly distributed 
and their spacing is much smaller than 2 J for pairs of levels. On the other hand, the spectrum does not 
show a spectral fragmentation, in contrast to the Bose- Hubbard model in a double well [33l |. where only 
in the high-energy part of the spectrum nearly degenerate pairs of levels appear. 

The difference between harmonic and anharmonic cavities is even more pronounced for the dynamics 
of the return and transition amplitudes. While there is only a periodic behavior with the single frequency 
J in Eq. (jlip . anharmonic cavities have a more dynamic behavior (cf. Figs. [5J|31). In particular, on the 
time scale considered in Figs. [51 [21 there is no periodic behavior but oscillations on much shorter scales 
than tt/ J sa 4. This is a consequence of the fact that the individual energy levels Ek — J{2k — N) in Eq. 
(jlOp are invisible in the dynamics of the harmonic cavities due to 

k=0 ^ ^ fc=0 ^ ^ 

Such kind of interference effect is accidental for harmonic cavities and does not occur for anharmonic 
cavities. Therefore, we can distinguish the individual levels in the dynamics only of the latter. 
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FIG. 1: Spectral density po,o of two coupled cavities with two-level atoms (coupling strength g « 1.2) for 100 
photons, inter-cavity tunneling rate J = 0.8. All energies are measured in units of the cavity frequency hujo. 



The return amplitude decays for both systems rapidly (cf. Fig. [2]) but it recovers much earlier for the 
anharmonic cavities, not to the full value though. Remarkable is the behavior of the transition amplitude. 
The time Tt it takes to reach the state |A^, 0) from |0, A^) for the first time is about the same for both 
systems, indicating that Tt ~ tt/ J must be solely determined by the tunneling rate J. This would allow 
us to measure the tunneling rate J in the dynamics of the system, regardless of the anharmonicity. 

Our main goal, the dynamical creation of an entangled state from a pure state, is also strongly af- 
fected by Tt, since entanglement in terms of a NOON state is not possible for times shorter than Tt- For 
times larger than Tt only the anharmonic cavities can reach the state \N, 0) while maintaining a non-zero 
overlap with the initial state (cf. Figs. [21 13]). On the other hand, only for a small number of photons 
(e.g., TV = 2) the harmonic cavities are capable to create a NOON state dynamically. For a discrete 
sequence of time steps we have counted the occurrence of certain values of the return and transition am- 



plitude ((iV,CT;0,cr| 



iHt 



N,a;0,a),{0,a;N,a\( 



iHt I 



N,a]0,a)) for harmonic and anharmonic cavities. 



This yields the conditional probability P(co,cjv), which was defined at the end of Sect. |lTl Example is 
plotted in Fig. 01 The plots demonstrate that the dynamical creation of a NOON state is feasible for 
harmonic cavities with up to iV = 4 photons while for anharmonic cavities this can be achieved even for 
N = 100. In comparison to bosons in a double well with Hubbard interaction ;33i] this probability is 
quite small though. 



V. DISCUSSION AND CONCLUSIONS 



We have considered a pair of optical cavities, where the photons in each cavity couple to a two-level 
atom. Then the cavities are described as JC models. For the initial state both cavities are prepared 
in an eigenstate of the JC Hamiltonian. Then we have connected the cavities by an optical fiber such 
that photons can tunnel between them. The resulting evolution of the quantum state of the combined 
system is determined by an effective Hamiltonian that resembles the Bose-Hubbard model with modified 
photon-photon interaction. The cavity frequency luq is assumed to be the same in both cavities. In this 
case Wo provides only a global shift of the spectrum, whereas the level spacing is entirely determined by 
the tunneling rate J of the optical fiber. For harmonic cavities (i.e. in the absence of the two-level atoms) 
the distribution of the levels with spectral weights pj = K^ol^-j)!^ is binomial with equidistant energy 
levels. The resulting evolution is periodic and corresponds to Rabi oscillations with a single frequency J. 
This behavior was observed experimentally for weakly interacting bosonic atoms [4], Q and should also 
be accessible for photons in harmonic cavities. The amplitudes for visiting the initial Fock state | A^, 0) or 
the complimentary Fock state |0, N) vary as cos^ (Jt) or (— i)^ sin^(Ji), respectively. This implies for a 




FIG. 2: Return amplitude as a function of time for two coupled harmonic cavities (full curve) and for two coupled 
cavities with two- level atoms (dashed curve). The parameters are the same as in the previous Figure, the time is 
measured in units of I/cjq. 
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FIG. 3: Transition amplitude as a function of time for two coupled harmonic cavities (full curve) and for two 
coupled cavities with two- level atoms (dashed curve). The parameters are the same as in the previous Figure. 



large number N of bosons that (i) these states are visited only for a very short period of time and (ii) the 
tvifo Fock states are visited at different times. Thus the dynamical creation of a NOON state from a Fock 
state |iV, 0) is very unlikely for harmonic cavities, unless the number of photons is small. The reason is 
that the photons can travel without seeing each other through the entire Hilbert space. A simultaneous 
overlap of \'i>t) with both Fock states \N, 0) and |0, N) is very unlikely then. This is a situation in which 
it is very difficult to control and follow the quantum evolution. On the other hand, applications of finite 
quantum systems, such as in quantum information processing |4ll. |42||. require a controllable evolution, 
in which only certain parts of the available Hilbert space can be visited with reasonable probability. 
In terms of our two-cavity system this means that the spectral weight pj = \{'^o\Ej)\'^ with respect to 
the initial state |^o) is small for most eigenstates \Ej) and has only a few pronounced maxima that 








FIG. 4: Probability P(co,cjv) for creating a NOON state in coupled harmonic cavities with A'^ — 2,4,6 (from 
right to left in the left panel) and in coupled anharmonic cavities with A^ = 100 (right panel). The parameters 
are the same as in the previous Figure. 



can be used for information storage. We have found that such a structured spectral density appears 
for anharmonic cavities, created by coupling two-level atoms to the cavity photons. Then the photons 
experience a mutual influence vifhich restricts their individual random walks in Hilbert space significantly 
and, what is even more important here, they can have a simultaneous overlap with both states |0,A^) 
and |0,A^). This effect enables the system to create dynamically a NOON state. The latter allows us 
to conclude that the complex quantum dynamics of two coupled anharmonic optical cavities offers an 
approach for quantum information processing as it has also been proposed for ultracold atoms [41| and 
cold trapped ions p3 |. 
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